Faceting and stress transfer in the missing-row reconstruction of Ir (110) 



OO 
On 
ON 



3 
^1" 



-a 
c 

o 
o 



> 
o> 

o 

OO 
On 

-i— > 
c3 



i 

C 

O 

o 



X 



Alessio Filippetti 1 ' 2 and Vincenzo Fiorentini 1,3 
^ Istituto Nazionale per la Fisica della Materia and Dipartimento di Fisica, Universita di Cagliari, 1-09124 Cagliari, Italy 

(2) Department of Physics, University of California, Davis, CA 95616, U.S.A. 
(3) Walter Schottky Institut, Technische Universitdt Miinchen, Am Coulombwall, D-85748 Garching bei Miinchen, Germany 

(February 1, 2008) 

We present ab initio total energy and stress calculations for the unreconstructed and (2x1)- 
missing-row reconstructed Ir (110) and Rh (110) surfaces. We use a model based on ab initio results 
to show that the (nxl) reconstruction is a faceting transition to a long-wavelength-corrugated (111)- 
like surface, and find the 3x1 structure to be the most stable for Ir in accordance with experiment. 
We then use the stress density to analyze the stress increase upon reconstruction, and find it to be 
due to a changed balance of tensile and compressive contributions in the near-surface region. 



Experimental findings! clearly support the missing-row 
reconstruction model for the (110) surfaces of the 5d met- 
als Ir, Pt and Au.B The reconstructed surfaces have been 
observed in a variety of n X 1 structures (presumably 
not all at equilibrium) whereby n-m out of n atomic 
rows of first-neighbor atoms, oriented in the [110] direc- 
tion, are removed from the m th plane, so that all planes 
through the (n+l) th expose rows to the vacuum. Calcu- 
lations support the existence of this reconstruction pe 
tern for Auou and Ira and its absence for Agu and Rh 
For Ir, mixed reconstructions patterns with predomi- 
nantly 3x1 character, have been observed,! while first- 
principle calculationsB confirmed that at least the com- 
putationally affordable 2x1 structure is indeed favored 
over the unreconstructed one. On a different line, swj-i 
face stress has baHjjnvestigated both experimentally^ 
and theoreticallya'BEj (mostly for Pt, Pd, Rh, and Ir). 
The behavior and role of stress in the reconstruction 
remains to be investigated. Also, the stability of nxl 
structures has not yet been supported theoretically from 
first-principle (accurate semiempirical studies exist on Au 
(110)0). 

Here we address the existence of the 2x1 reconstruc- 
tion by means of ab initio calculations. We then set up 
a model to describe nxl reconstructions, based on ab 
initio quantities for the computationally tractable cases. 
Finally, we analyze the behavior of surfacej-siress using a 
recently-developed tool, the stress density.til In the first- 
principles local-density-functional calculations we use ei- 
ther fully separable norm-conserving (Ir) or ultrasoft 
(Rh) pseudopotentials, a plane-wave basis cut off at 40 
Ry (Ir) and 30 Ry (Rh), and surface slabs of up to 11 
atomic planes, fully relaxing all structures. For further 
details, see Ref. g. 

Ab iaitio results -.— .In Table || we report surface 



t iHTgylLHl surface stressEl, and workfunction of the unre- 
constructed Ir and Rh (110) surfaces, and their changes 
upon reconstruction. The reconstruction energy H lcc 
is the difference between reconstructed and unrecon- 
structed surface energy, i.e. H lcc < means favored 
reconstruction. This is found to be the case for Ir, not 
for Rh, in agreement with experiment nu For Ir we find an 
energy gain of 0.03 eV/atom, close to the ~ 0.05 eV/atom 



obtained 



previousl;ylill 



for Au. 



The unreconstructed surfaces are under a large tensile 
stress (negative in our convention), as already .custom- 
ary from previous work on transition metalsJj'EMLj The 
stress is appreciably larger fopJr than for Rh due to the 
influence of relativistic cffectsEiJ) , and highly anisotropic: 
A=T u nre cy r unrec i s 1.89 for Ir and 1.61 for Rh (here x = 

[001] and y = [110]). The anisotropy is almost entirely 
due to relaxations (in the ideal structure we find A=1.17 
and 1.11 for Ir and Rh), in agreement with Feibelman's 
argumenta that the first-neighbor bonds having a com- 
ponent along the [001] direction can shorten upon re- 
laxation, and relieve a large part of their tensile stress, 
whereas first- neighbor bonds forming the [110]-oriented 
atomic rows cannot shrink because of symmetry con- 
straints, and thus conserve the stored tensile stress. Our 
A-values are similar to Feibelman'sEI for Pt and Pd. Both 
A and its change upon relaxation are greater for the 5d 
elements - that is, the larger the surface stress, the larger 
the portion of stress relieved by the relaxation. 

Reconstruction causes a puzzling behavior of the sur- 
face stress and workfunction. As a tensile stress is nega- 
tive by convention, negative Ar's such as those in Table || 
indicates an increase in tensile stress. The puzzling result 
is thus that the surface stress is more tensile on the recon- 
structed surface than on the unreconstructed one. This 
feature is surprising, as the surface stress is mostly stored 
into the bonds of the [110]-oriented rows, and one would 
expect it to decrease upon removal of one row out of two. 
In addition, the workfunction increases with reconstruc- 
tion, suggesting a decrease in surface roughness. This is 
surprising because reconstruction decreases the surface 
atomic density and apparently increases the roughness 
(the corrugation amplitude is doubled). (Although this 
is not our concern here, we note that stress does not de- 
crease upon the.-2 x 1 transition, which seems to disagree 
with arguments!^ to the effect that stress is the driving 
force of this transformation. Stress has been invoked to 
explain e.g. the mesoscopic patterning of period ~ 1500 
A recently observed!] on 2x1 Pt (110): that, however, 
is a further transformation of the 2 x 1-reconstructed sur- 
face.) 
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The faceting model - The increase in workfunction 
and surface stress with reconstruction can be explained 
easily if we view the reconstruction as a microfaceting 
procesp,|iato (lll)-like facets to vacuum. In previous 
works,Eni-3 faceting had already been suggested as a key 
feature of the missing-row reconstruction: using our first- 
principle results, we now set up a model giving quanti- 
tative evidence in favor of the faceting hypotesis. We 
then extend it to arbitrary n. The reconstructed surface, 
sketched in Figure |l|, has (lll)-oriented facets (shaded 
areas) at an angle of ~ ±35° from the (110) plane. The 
average coordination number and bonding geometry of 
the facet atoms are the same as on the (111) face, sug- 
gesting that the facets may exhibit properties close to 
those of the (111) surface. 




FIG. 1. Faceting model geometry for the missing-row re- 
construction: the shaded area is the quasi-(lll) reconstructed 
surface. 

We verify such a similarity in a simple but stringent 
way, comparing directly-calculated first-principles ener- 
gies and stresses for the reconstructed surface with a 
model assuming the energy and the stress per unit facet 
area to be equal to those of the relaxed (111) surface 
(also calculated ab initio), and integrating these values 
over the facet area: 
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Here do is the bulk lattice constant, d the first-to-third 
plane distance (see Figure |l|), cr(m) M(m) tne ener gy 
per unit area of the relaxed (111) surface. Analogously, 
for the surface stress we take, 
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and for r™ c the projection of t^ii) along x. For the 

unrelaxed surface, d — ao/\/2, and er rcc = 2 07m). In 
Table || we compare the ab initio and model-derived val- 
ues for the relaxed surfaces. The agreement of the two 
independent evaluations is indeed very good. It is then 
plausible to identify the 2x1 reconstructed (110) surface 
with a spatially modulated 'quasi- (111)', quite similar to 
the true (111) surface. The effect of outward and inward 
edges appears to nearly compensate each other on aver- 
age, and the reconstruction can be thought of as a transi- 
tion from a rough surface (the unreconstructed (110)) to 



a close-packed, high-coordination surface, i.e. the quasi- 
(111) surface, perturbed by a relatively long-wavelength 
corrugation (A=2ao). Thus, the increase of workfunction 
is justified as it actually reflects a transition from a (110)- 
like to a (lll)-like situation (the behavior of the surface 
stress is discussed later). 

The model easily explain the different behavior of Ir 
and Rh: the energetics of reconstruction is governed by 
the ratio R — cr(uo) /< T (m) • From Equation (|l|) we have 
that the missing-row reconstruction occurs if 
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where ^4(m) = v3oo/4. Therefore, the surface recon- 
structs if R is larger than a reconstruction threshold Rt 
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where Ad is the change of d from its ideal value d ld = 
ao I V%, and the second equality is valid to first order in 
(Ad/d). Rt is in other words the ratio of the facet area 
(a function of d) to the (111) area. In Table 111 we re- 
port values of R, Rt, and Ad/d lA for Ir and Rh. If we 
disregard relaxations (i.e. put Ad = 0), then Rt=2, and 
neither Ir nor Rh should reconstruct. The reduction of 
d (and hence of the facet area) is the essential ingredient 
for Ir reconstruction. Rh, on the other hand, has nearly 
equal Rt, but a much smaller R, i.e. a much smaller 
surface energy anisotropy. Notice however that Rh is ac- 
tually just on the verge of reconstruction; according to 
the model, Rh would also reconstruct if it had the same 
relaxations as Ir (i.e. if the threshold values for Rh and 
Ir were the same). 

An n-dependent faceting model - Th model just dis- 
cussed is n- independent. Since 3x1 faceted domains 
have been preferentially observed in experiment Jj we 
need to extend the model to a general n x 1 recon- 
struction. The vertical and horizontal distances between 
edges and valleys are now given by djf cal = na /2y2 
and dhoriz = nao/2. The surface energy (per lxl area) 
generalizes to 
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where d n is the distance between first and (n+l) th layer; 
the threshold ratio becomes 
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with Y^,a Pi = Sd n /d lAeli[ the sum of relaxations expressed 
as fractions of the ideal interplanar distance. The recon- 
struction is favored for 
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Since the p^s are not known for layers far away from 
the top, we assume hypotetical relaxations extrapolating 
to zero atpjthe tenth plane below the surface the values 
calculated! 5 ! for the 2x1 case. Ia^m) is displayed as cir- 
cles in Figure @: the reconstruction is favored for n < 4. 
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2. The (110)-n x 1 reconstruction energy predicted 
faceting model as a function of n (see text). 



Explaining the stress increase: application of the stress 
density - In the light of the faceting model, the ratio- 
nale for the increase in tensile stress is that reconstruc- 
tion, while removing rows from the surface, exposes new 
ones to the vacuum, causing the stress stored in the sur- 
face rows atoms to be spatially redistributed, and specif- 
ically, transferred towards the substrate. To support this 
argument directly, we introducedjld in analogy to the 
Chetty-Martin energy densityjia the stress density, de- 
fined as a tensor T a p(r), functional of the charge density, 
whose integral over the cell equals the macroscopic stress: 



drT a p{r) = t, 



'a/3 ■ 
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The detailed formulation, as well as a discussian of gauge 
dependence issues, will be given elsewhere.til The core 
idea is that while surface stress only provides informa- 
tion on certain phenomena at surfaces in an integrated 
fashion, the stress density allows to analyze locally its 
different components. In the present case, we need to 
consider the planar average T a p{z) of the stress density 
along the [110] direction, and the average of T a p{z\ over 
the inter-layer distance (the macroscopic averagetj). In 
Figure |^, we compare planar and macroscopic average of 
the stress density of 



The monotonic decrease of ioVm) as n— > 1 disagrees with 
two facts: a) the lxl should be unstable with respect to 
the 2x1 (although one may simply invoke the inapplica- 
bility of the model for n—1); b) experiments^ show that 
the 3x 1 is more stable than the 2x1. To resolve this, one 
has to account for the contribution of the edges. Adopt, 
ing a model due to Marchenkota and Shchukin et alx3 
based on elasticity theory, we add an edge energy and an 
elastic energy to the reconstruction energy I a (my The 
total reconstruction energy is given by 

C\ C<z 
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The parameter C\ is taken to be equal to la/ my he. 
the energy needed to create on a (111) face the edges 
and valleys of a (HO)-nxl. For C2, we useE3 



C7 2 = 8(1 - v 2 ) 
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with f/>=0.611 rad ~ 35°|— the faceting angle, Y—5.2J& 
Mbar the Young modulus,l!3 v — 0.26 the Poisson ratio,tZl 
and ao=7.289 bohr the calculated lattice constant of Ir, 
T (ni) — 1-96 eV/ (lxl) the first-principles calculated sur- 
face stress of Ir (111). The successive approximations are 
depicted in Figure g. The circles are io'(iii); adding the 
edge energy we obtain (squares) a minimum at n=2, and 
stable n=3 and n=4; the complete expression (filled di- 
amonds) gives n=2,3 stable with n=3 lowest in energy. 
The calculated reconstruction energy is not reproduced 
quantitatively, but its qualitative behavior is now correct. 
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FIG. 3. Planar and macroscopic averages of stress density 
for the reconstructed (full lines) and unreconstructed (dashed 
lines) Ir (110) surface. 

the unreconstructed (dashed lines) and reconstructed 
(full lines) surfaces. Both surfaces are considered in their 
unrclaxcd configuration for graphical convenience and, 
exploiting inversion symmetry, we only display one slab 
half. Note that the two inequivalent stress components 
along x and y are indistinguishable on this scale. 

The macroscopic averages vanish into bulk (by the 
equilibrium condition of vanishing stress) and in the vac- 
uum, and they deviate from zero whenever surface effects 
arise. Thus, they enable one to establish the effective 
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boundaries of the surface region as far as stress propa- 
gation is concerned. The planar averages show periodic 
oscillations in the bulk region, and vanish into the vac- 
uum. The negative peaks are tensile contributions arising 
mostly from electron-ion interactions, and are located at 
the atomic positions. The positive peaks rising between 
atomic planes store the compressive stress of essentialy 
kinetic origin. 

The stress transfer upon reconstruction is easily de- 
tected moving from left to right in Figure ||: for 
the unreconstructed surface (dashed line), a non-zero 
macroscopically-averaged stress density is only visible in 
the first layer, whereas for the reconstructed one it starts 
off to non-zero values one layer earlier. The stress is 
indeed redistributed from the surface to the substrate. 
Focusing on the planar averages, the negative (tensile) 
peak in the first layer of the reconstructed surface is re- 
duced with respect to the unreconstructed one. This cor- 
responds to the expected reduction of tensile stress upon 
removal of the [110] -oriented atomic rows. On the other 
hand, the positive peak between surface and substrate is 
also strongly reduced, and it compensates (in fact, over- 
compensates) the reduction of the tensile component. In 
other words, row removal reduces the tensile stress in the 
surface plane, but at the same time it allows electrons 
to delocalize into the troughs thus created- reducing the 
compressive kinetic stress. Ho and BohenQ arrive at es- 
sentially the same conclusion by analyzing the individual 
contributions to the total energy, and finding a reduction 
in electron kinetic energy, larger in Au than in Ag. 

In summary, we showed that the (110) missing-row re- 
construction can be understood as microfaceting of the 
(110) into a bent, quasi-(lll) surface. The reconstruc- 
tion reduces the effective surface roughness and causes 
a stress transfer from the surface towards the substrate. 
The latter effect has been analyzed by means of the stress 
density. The different behavior between 5d and 4d ele- 
ments can be traced back to the greater (110) /(111) sur- 
face energy anisotropy, and ultimately to the relativistic 
character of the heavier elements.L3 
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TABLE |. Surface energy a unrcc , stresses r^ nroc and 
Ty nrcc , and workfunction w unrcc of unreconstructed 
(110); H rec , AW, At x and Ar y are changes in previ- 
ous quantities upon reconstruction. Results in eV or 
eV/atom. 



^.unrec 


H ICC W unrcc AW T-r rcc At* r H unroc Ar y 


Ir 2.59 


-0.03 5.42 +0.15 -1.70 -0.44 -3.21 -0.50 


Rh 1.89 


0.07 5.07 +0.15 -1.25 -0.09 -2.01 -0.13 


TABLE [El]. Reconstruction energies H rcc and recon- 
structed surface stress r r (eV/atom) as calculated di- 
rectly (c) and by the faceting model (m) (see text). 


Ir Rh 




surt surt rrrcc surt surt 
' x 'y 'x >y 


c -0.03 


-2.14 -3.71 0.07 -1.34 -2.14 


m -0.04 


-2.20 -3.82 0.03 -1.32 -2.28 


TABLE [II 
relaxation c 


■ R ~ <T [iio]/< T [iii]i faceting threshold Rt', 
hange of first-third layer distance Ad/d ld . 


R R T Ad/d id 


Ir 


1.977 1.947 -0.078 


Rh 


1.948 1.953 -0.069 
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